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Abstract 

We present the results of a numerical calculation of the photon splitting 
rate below the electron-pair creation threshold {lo < 2m) in magnetic fields 
B ^ Bcr = rr? je = 4.414 x 10^ T. Our results confirm asymptotic approxima- 
tions derived in the low- field {B < Bcr) and high-field {B ^ Bcr) limit, and 
allow interpolating between the two asymptotic regions. Our expression for 
the photon splitting rate is a simplified version of a formula given by Mentzel 
et al. We also point out that, although the analytical formula is correct, the 
splitting rates calculated there are wrong due to an error in the numerical 
calculations. 
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The exotic process of magnetic photon sphtting, i.e. the decay of a photon into two 
photons in the presence of a very strong magnetic field, has recently attracted renewed 
attention, mainly because of the great importance this process may have in the interpretation 
of the spectra of cosmic 7-ray burst sources. The basic formulae for magnetic photon 
splitting had already been derived in the seventies In the first approach the photon 

splitting effect was analyzed using the Heisenberg-Euler effective Lagrangian. This method 
is justified under the condition u; <^ m. The result was the {uj /mY{B / Bcr)^ dependency 
of the attenuation coefficient in the weak field regime {B <C B^r) IH],!!- Adler was the 
first to solve the problem for arbitrary magnetic field strengths and photon energies below 
the pair creation threshold. He used the gauge invariant proper-time method and presented 
numerical results up to -B = B^r for the two cases uo = m and uj ^m. Other gauge-invariant 
versions of the splitting-amplitude were found later by Stoneham 0, Baier et al. 0, and 
recently, using a path- integral approach, by Adler and Schubert 0. Nevertheless, numerical 
results, apart from those obtained by Adler, and in particular for magnetic fields exceeding 
Bcri were not available. Mentzel et al. |§] therefore undertook a rederivation of the photon 
splitting rate using a configuration space representation for the electron propagator in a 
strong magnetic field with the aim to find an exact analytical expression that would be 
suitable for numerical evaluation for magnetic fields above B^r- The representation chosen 
for the electron propagator, though not gauge- invariant, had the advantage that the final 
expression for the photon splitting rate did not contain any singularities for photon energies 
below uo = 2m. Moreover, the sums over the Landau-states that had to be calculated in this 
method were expected to converge very rapidly for B > B^r- 

The numerical results presented in seemed to suggest that the absolute values of the 
rates were much bigger than was previously assumed, but failed to reproduce the asymptotic 
low-frequency, low-field behavior that followed from the earlier investigatons. Consequently, 
these results were criticized by Adler , who suggested that the fixing of the gauge had led 
to calculation errors that could not be detected in the final expression. 



More recently, Baier et al. [10| determined the behavior of the photon splitting rate in 



the limit of extremely high magnetic fields, B / B^ ^ 1, and demonstrated that in this case 
the splitting rate becomes independent of field strength. In agreement with these results 
are the recent results of Heyl and Hernquist [|ll|] who reexamined the effective Lagrangian 
method and were able to calculate the attenuation coefficient for arbitrary magnetic field 
strengths and u <^ m. 

In this note we will reexamine the formula given by Mentzel et al. simplify it, eval- 
uate it numerically, and compare with the results of the asymptotic formulae obtained by 
Adler P[, Baier et al. |]TU|, and Heyl and Hernquist |TT|. This will allow to quantitatively 



assess the range of validity of the different asymptotic results. We point out that the nu- 
merical results of [H are in error because of an incorrect sign in the code for evaluating the 
analytical formula for magnetic photon splitting. 

Following the paper by Baier et al. we shall concentrate on the process -L^|| ||, where 
a photon with the polarization vector parallel to the plane formed by the k vector and the 
magnetic field axis decays into two photons with polarization vectors perpendicular to that 
plane. This is expected to be the dominating decay channel as soon as dispersive effects 
become important. Note that these designations for the polarizations are exactly opposite 
to those chosen by Adler. Furthermore, we confine ourselves to photon energies below the 
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pair creation threshold uj < 2m, above which pair creation is expected to dominate photon 
sphtting. 

The complexity of the general expressions for the photon splitting rate in magnetic fields 
has always been an impediment to quantitative studies of the importance of the photon 
splitting process in actual astrophysical model calculations. The availability of appropriate 
approximation formulae, and a knowledge of their respective ranges of applicability, is there- 
fore a prerequisite for such investigations |T^. On the low-frequency, weak-field {B ^ Bcr) 



side, the attenuation coefficient can be written (in the notation of Stoneham 

„ 1 , MUX a^m / B \ ^ / Lu 



where the scattering amplitude (Mi(_B)) is given by eq. 41 of ref. p. Adler |^ demonstrated 
that for field strenghts up to ~ 0.1 Bcr the coefficient is excellently reproduced by that 
equation, with Mi replaced with its = value, Mi(0) = 26/315. 

On the very high field side, B ^ B^r Baier et al. |Ty] have recently shown that, for 
given energy below 2m, the attenuation coefficent tends towards a constant, magnetic field 
independent value, viz. 

do;' 



^-B>Bcrl (2) 



327r' ^2 

where Tb-^Bc^ is the high- field approximation of the splitting amplitude: 



2(47ra) 



7b>b„ — Y 



: arctan 



uj"\/Am? - uj""^ \ ^/Arrfi - uj""^ 



H = arctan — uj . (3) 

cuV4m2 - VV4m2 J 4 J ^ ^ 

On the basis of an accurate numerical evaluation of the general expression for the photon 
splitting rate we will now determine the range of field strengths and frequencies from where 
on the high-field asymptotic behavior is valid, and provide data for the range of intermediate 
field strengths. 

The general expression we evaluate is the one derived in [Q. Two important remarks on 
this work are in order. The first is that by reason of a simple error in sign in the computer 
code (the quantity K2 defined in eq. (A12) in the Appendix of ref. contained a plus sign 
in the code instead of a minus sign in the first term on the right-hand side) all the numerical 
results are wrong. This error escaped the numerical checks of the polarization selection rules. 
Speculations that gauge invariance problems caused the deviations from previous results are 
thus found to be insubstantial. The second remark is that the analytical expression obtained 
by Mentzel et al. P] is correct. The low-frequency, weak-field limit (||) can be derived from 
it analytically, and the results of the numerical evaluation (corrected for the error in sign) 
are in agreement with previous results in all ranges of parameters where a comparison is 
possible. We note that by a suitable rearrangement of terms the numerical evaluation of the 
expression of could be speeded up by two orders of magnitudes. The main step in these 
optimizations is the summing up of the different polarization states. We performed this sum 
using the computer algebra program Maple. Once the summation has been carried out, the 
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remaining integration can be reduced to 5 integrals that can be evaluated independendly 
from the rest of the formula. 

An important point with regard to the accuracy of the results is the inclusion of suf- 
ficiently many Landau states in carrying out the sums over the intermediate states. In 
particular for magnetic fields in excess of Bcr these sums are found to converge very rapidly. 

In Fig. 1 we present, for the three frequencies uj/m = 0.1, 1.0, 1.99, the results of the 
numerical evaluation (crosses) of the general expression for the magnetic photon splitting 
attenuation coefficient for the decay channel _L— >■ || || as a function of the field strength in 
comparison with the predictions of the asymptotic formulae. The steep solid line on the left 
corresponds to the dependence of the low-field formula eq. (|1]) (with Mi(0) inserted for 

Mi(i?)), while the dotted horizontal line corresponds to the high-field limit @). The figure 
shows how, for all three frequencies, the numerical results approach the low-field limit as B 
decreases, and that the numerical values tend towards the high-field limit as B increases. 
The results are in good agreement with those of Heyl and Hernquist |TI]. We find that their 
formulae have a very high accuracy for photon energies up to = 0.5m. The high-field limit 
of Baier et al. is reached, to within a few per cent, already for magnetic field strengths 
B / Bcr around 30, independent of the photon-energy. 

We now proceed to results for the differential splitting rate. This quantity provides 
the information on the probability with which a certain partitioning of the energy of the 
decaying photon on the outgoing photons is assumed. Therefore it is reasonable to normalize 
this quantity to unity. Note that in this respect our procedure differs from that of Baier 
et al. |T^, who used a magnetic-field dependent normalization constant, which blurs the 



comparison of the results at different field strengths. 

In Fig. 2, for the energies uj/m = 0.1, u/m = 1.5 and u/m = 1.99 and magnetic 
field strengths B = B^r, 5 B^r and 100 B^r numerical results are shown for the differential 
splitting rate (normalized to unity when integrated over the abscissa) as a function of uj'/uj. 
For B = 100 Bcr, as expected, the shape of the differential splitting rate is exactly identical, 
for all three photon energies, to the one determined from Baler's asymptotic expression 
(solid curves in Fig. 2). This agreement, however, is seen to be present even for the smaller 
magnetic field strengths of i? = 5 B^r and B = B^r, and it is only for photon energies 
approaching the pair creation threshold in the latter case that sizeable deviations from 
the asymptotic behavior appear. We also see from Fig. 2 that the plateau-forming effect 
described by Baier et al. |10| occurs already below B = 5 B^r- Thus we have the result that 
the high-field asymptotic formulae represent a good approximation even at field strengths 
where they were not a priori expected to apply. This opens the possibility of carrying out 
astrophysical model calculations using the asymptotic high- field formulae for B ^ bB^r, field 
strengths that are supposed to be present in soft 7 repeaters [|l^,|l^]. 

In this note we have presented numerically accurate values of attenuation coefficients 
for photon splitting in strong magnetic fields, and have used these results to determine the 
ranges of parameters in which asymptotic approximation formulae can be applied in the 
place of the full complicated expression for the photon splitting rate. The results should 
stimulate, in high-energy astrophysics, new quantitative studies of the role of the exotic 
process of photon splitting in the formation of high-energy spectra of strongly magnetized 
cosmic objects. 
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FIG. 1. Total attenuation coefficient for tfie splitting channel _L^|| || as a function of field 
strength for the three photon frequencies u/m = 0.1, 1.0, and 1.99. Crosses denote the results of 
the present calculations, the weak-field result is represented by the solid straight line, the high-field 
approximation by the horizontal dotted line. 
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FIG. 2. Differential splitting rate (normalized to unity when integrated over the abscissa) for 
the splitting channel \\ and three frequencies u of the incoming photon as a function of the 

ratio of the frequency uj' of one of the outgoing photons and uj, for the magnetic field strengths (from 
top to bottom) B/Bcr = 100, 5, and 1. The solid curves represent the results of the asymptotic 
high-field approximation by Baier et al. This approximation is seen to be valid even at magnetic 
field strengths as low as S ~ 5Bcr. 
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